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Abstract 

In this paper we present a multilayer particle deposition model on a random tree. 
We derive the time dependent densities of the first and second layer analytically and 
show that in all trees the limiting density of the first layer exceeds the density in the 
second layer. We also provide a procedure to calculate higher layer densities and prove 
that random trees have a higher limiting density in the first layer than regular trees. 
Finally, we compare densities between the first and second layer and between regular 
and random trees. 

AMS 2000 subject classification: 82C22, 82C23. 

Key-Words: Car parking problem, Random sequential adsorption, Sequential fre- 
quency assignment process, Particle systems. 

1 Introduction 

Parking models with screening were first studied in the field of ballistic particle deposition, 
see for example [I]. In those models particles are moving towards a substrate or a fiber until 
they encounter a previously deposited particle or the substrate itself. A particle always 
tries to park on a layer as low as possible but due to "screening" the particle cannot pass 
formerly deposited particles. In our model the screening rule makes every particle park in 
the highest layer possible where it is supported by a particle in the layer below (see figure 
1). This "Tetris" model is very different from the so-called Sequential frequency assignment 
process (SFAP) in which particles (assignments) can skip particles on their way down [2]. 
In that model particles are deposited in the lowest layer (frequency) possible. It has been 
found for the SFAP that there is an increasing limiting density of particles in higher layers 
due to boundary and other effects [SIE]. In this paper we will show analytically that in 
the model with screening the opposite is true. The density in the first layer appears to be 
higher than in the second layer. We conjecture that the same applies to the other layers. 
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Also, we generalize the model on the one- dimensional lattice to a model on regular and 
random trees. Recently, several random particle deposition tree models have been studied 
in [4JEJE]. However, to our knowledge this is the first time a multilayer random tree model 
is treated. 



Multilayer Parking with Screening on a Line 
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Figure 1: Example of a realization of the particle deposition process with screening on a 
one-dimensional lattice. 



2 Layer Densities 
2.1 The Dynamics 

The precise definition of the model is as follows. We consider a random tree with vertices 
i G V and degree at the site i given by Di. We choose Di to be independent random 
variables with the same distribution Q given by 

Q(A = k) = a k (2.1) 

on the integers starting from 2. The latter requirement ensures that we have no open ends 
with probability one. 

We denote the generating function of the distribution by 

oo 

G{s) = Y J a kS k (2.2) 

k=2 
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We will denote the expected value with respect to this probability distribution by the 
same symbol Q. We fix a realization of the random tree, and, denoting by V its vertex set. 

In the first part of the paper we will be interested in the behavior of particle config- 
urations arising from particle deposition in the first two layers, that is we consider the 
marginal of an infinite layer particle model on the first two layers. This will be generalized 
to higher layers later. To describe the behavior on the first two layers we consider (suit- 
ably coded) occupation numbers m = (m(i))i e v £ ^ = {0, 1,2,3} V . Here the spin m(i) 
denotes the joint occupation numbers at vertex i at height 1 and 2. It is useful for short 
notation to interpret the occupation numbers at various heights as binary digits and write 
ordinary natural numbers. That is we write 

{0 if vertex i is vacant in the first and second line 

1 if vertex i is occupied in the first but not in the second line , . 

2 if vertex % is occupied in the second but not in the first line 

3 if vertex i is occupied in the first and in the second line 

so that m(i) G {0,1,2,3}. 

Now we can describe the dynamics of the joint process of particle occupations in the 
first two layers and the total number of particles which have arrived by the following 
generator. Note that it is really necessary to consider also particle arrivals beyond the 
first two layers because of the screening effects higher layer particles might have on lower 
layers. It is furthermore necessary to distinguish two sorts of particle arrivals: those which 
change the lower layers and those which leave the lower layers unchanged. Let F be a 
joint function of particle occupations in the first two layers and particle numbers. Then 
the generator of our process reads 



£F(m,N) = J2( E F(rn s >\N k )r k (s;<m k ) 

keV s=l,2,3,4 

+ (1- Yl r k (s;Wl k ))F(m,N k )-F(m,N)^ 



(2.4) 



with 



and with 



[ m(i) it k 7^ i 

where 9Jl k := (N(l), m(l))ie{k}uC(k) where C(k) := {i : dist(i,k) = 1} is the neighborhood 
of vertex k. 

So, the first term in the generator describes the events when the addition of a new 
particle also changes the configuration in one of the first two layers. The second term in 
the generator describes the events when the first two layers are already full or screened, 
and a further adding of a particle does not change its filling. 
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The rates are either equal to zero or one. They are 1 precisely in the following cases 
(see also Figure [2]) listed below. 

1. i — > 1 Adding a particle in the first line at vertex i. We have 

r(l; {N(i) = 0} n {V ieC(i) : N(j) = 0}) = 1 (2.7) 

Indeed, this occurs when the site and all its neighbors are empty in all layers. 

2. i — >• 2 Adding a particle in the second line at i while the first line was empty at the 
site 



r(2; {N(i) = 0, 3 J C C(i) : V j E J : m(j) = 1, N(j) = 1, VJfe € C(i)\J : 

N(k) = 0}) = 1 



(2.t 



It is only possible to reach the state m(i) = 2 when at least one of the neighbors has 
a particle at layer 1 and all the others are totally empty. 

3. 1 i — > 3 Adding a particle in the second line while the first line was full at the site 

r(3; {m(i) = 1, N(i) = 1} n {Vj- eC(i) : JVtf) = 0}) = 1 (2.9) 

To get into state m(i) = 3 there must be one particle in vertex i and all neighboring 
sites should be empty to avoid screening. 

All other transitions are impossible. 
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Figure 2: Neighborhood configurations of a vertex that allow the mt transitions — > 1, 
l—s-2 and 1 — > 3 respectively. In this example the central vertex has three neighbors. The 
states are denoted with the notation (Nt,mt), e.g. (1, 1) means that one particle arrived 
and that it was deposited on the first layer. In order to have a transition from to 1 every 
vertex in the neighborhood has to be totally empty (Nt = 0). To get a transition from 
to 2 the vertex itself must be empty, but at least one of the neighbors has to have one 
particle in total that lies on the first layer. Finally, to have a transition 1 to 3 the vertex 
has to have exactly one particle located on the first layer, while the neighbors should be 
empty. 
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This generator defines a time-homogeneous Markov jump process on the infinite graph 
by standard theory [7] such that JZU) f t \ t =o^ m,N F(m t , N t ) = CF(m,N). 

Here K m ' N denotes the expected value with respect to the process, started in the initial 
configuration (to, iV) = (m(i), N(i))i £ v at t = 0. 

We underline that we consider the marginal of the first two layers of a model where 
particles may pile up to arbitrary high layers. The present model differs from the model 
discussed in [3j where particles that cannot be deposited in the first or second layer are 
rejected. 

2.2 Regular Trees 

We first consider the densities, taken at an arbitrary vertex called 0, 

pf(l) = P t (m(0) = 1) + P t (m(0) = 3) 
p d t {2) = P(m(0) = 2) + P t (m(0) = 3) 



on the first and second layer on a regular tree with degree d > 2, Having understood their 
behavior on a regular tree we can derive the densities on random trees easily in section I2T31 

Theorem 1 Consider the regular tree T d with degree d > 2. Particles arrive at the vertices 
of Td according to a Poisson process and obey the screening rules of deposition. Then the 
time dependent densities are, on the first layer 

1 _ e ~(d+l)t 

"< (1 > = 5+1 <"« 

and on the second layer 

^ - ( d \ d S-( d \ ^{-l)* d , d .-frW 



ft(2) \d-l) f^ Q \k) {d-l)k + d+l (d+1) 2 + (d+l) 26 

— V V s ( d \ d ~ fc (~ 1 ) fc r -\(d-i)k+d+\]t __+ p -(d+i)t 



e -[(a-i)K+a+x\t _ tc -(d+pt ^ 12) 

d-lj ^ \k) (d- l)k + d+ 1 d+1 



Proof of Theorem [TJ 

Throughout the paper we use the notation 

Df(s) = p(m(0) = s) 

for all s. We fix at a certain vertex 0. The surrounding vertices are numbered 1,2, ... ,d. 
First we calculate the time derivative of Df(T) and integrate back. Taking into account 
the first and third process depicted in Fig. 2 we see that 

Df(l) = P t (N t (0) = Q,V k < d N t (k) = 0) - P t (N t (0) = l,V k < d N t (k) = 0) 
= e ~(d+i)t _ te -t e -dt 

= -(t-l)e-( rf+1 )* (2.13) 
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So, now with D#(l) = we find 



D d (l) - - - e ~ (d+1 ^ + 1 tc~ (d+1)t 

t[) ~ {d + iy (d + 1) 2 + d + l 



(2.14) 



We apply the same technique to find Df(2) 

d 



D d t (2) = £ (^) P ^ Nt ^ = °^i<k N t(i) = hmtii) = ^k< 3 < d N t (j) = 0) 
d , ,x 

= Y ly k )Pt(Vi< k N t (i) = l,m t (i) = l,V k<j < d N t (j) = 0\N t (0) = 0)e~ t 

= e-*^Q[P t (iV t (l) = l,m t (l) = lliV t (0) = 0)] fe [P t (iV t (0) = 0)] rf - fc (2.15) 

Now we need to calculate the quantity Sf := Pf(N t (l) = l,m t (l) = l|iV t (0) = 0). This is 
done by constructing another differential equation. 



Sf = P t (N t (l) = 0,V i:d{iA)=1 N t (i) = 0\N t (0) = 0)-S? 



= e~ dt -S? 



whose solution is 



So, we have 



1 _ e -(d-l)t 
1^1 



(2.16) 



(2.17) 



Df(2) 



£ 

k=l 



-(d+l)t 



E 

k=l 



-(d+l)t 



1 + 



J_ e -t(l_ e -(d-i)t) 



1 _ e -(d-l)< 
1^1 



1 _ e ~(d-l)t 



-tid-fc 



d-1 



= _ e -(d+l)t + g -(d+l)t 



1 



-(d-i)t 



d-1 d-1 

C( / ,N / 1 N d-fc 



. _ e -,«». + e -, +1 ».gQ(_4_)- (_L) 

' h — n N / 



,-(d-l)kt 



3 -(d+l)t 



fcg-[(d-l)fc+d+l]t 



(2.18) 
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which gives (with £$(2) = 0) 



k) (d- l)k + d+ 1 d + 1 



I 1 c -(d+Dt _ ( d \ ( d \ d [{d - 1)k+d+1]t (2 lq) 

+ d+l 6 \d-lj f^ Q \k) (d-l)k + d+i e l2 ' iyj 

Finally, for Df(3) we find 

Df(3) = P t (N t (0) = l,V k < d N t (k)=0) 

= te- {d+1)t (2.20) 

so, that (with D$(3) = 0) 

Vt[6) ~ (d + i)2 (d + \y e d + i e [2 - 21) 

The densities of the first and second layer follow immediately by adding D d (1) and D d (3) 
for the first layer, and Df(2) and Df (3) for the second layer. □ 

Remark: Note that for the derivation of the formula for the first layer we did not have to 
use the absence of loops in a tree. Therefore, the first layer density on a graph is the same 
as on a tree, no matter whether they are regular or random. 



2.3 Random Trees 

Let us now consider the case of particle deposition on a random tree where the number of 
neighbors of every vertex is a random number according to some G(s) = X]^L 2 a n s" ■ We 
now have the following 

Theorem 2 Consider a multilayer random tree Td with generating function 
Gt(s) = E^^L 2 a ra s n . Particles arrive at the vertices ofTjj according to a Poisson process 
and obey the screening rules of deposition. Then the tree-averaged time dependent densities 
are, on the first layer 

00 (l- e -(k+l)t\ 

Opt(i) = Y. a *- — rxi — ( 2 - 22 ) 

k=2 

and on the second layer 

+ 1 «* E (t) E (*) (-»' (f ^t) I to**** (-4) 

d =2 fe=l V 7 i=0 V 7 \d=2 / JU 

where Zt := YZ2*de- {d -^/{d - 1). 
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Proof: 

First, we calculate QD t (l) and QD t (3). Notice that the derivatives of these functions in 
a certain vertex are not affected by the tree ensemble beyond the nearest neighbors. 
Therefore, we can immediately start averaging A(l) and A(3) over Q rather than dealing 
with its derivatives first. In the previous section we already found 

n d (l) - d d c -(d+i)t , 1 tc -(d+i)t (22^ 

Ut[l) ~ (d + i)2 (d + iy e + d + i [2 - 2b) 

where d now denotes the (random) number of nearest neighbors of the site under consid- 
eration. Averaging over Q results then in 

Similarly, we find 

OO OO / \ 



Adding these two results gives the density on the first layer. 
In the previous section we already found 

A(2) = e"*X; [P t dfc = l,m t (l) = l|JV t (0) = 0)] k [P t (N t (0) = 0)] d «~ k (2.28) 

where the number of neighbors of vertex i is denoted by di. Note that in this section the 
di's may be different since we are treating a random tree. So, we get 



do /j \ 

A(2) = e-'2 I )[S? k ] k [P t (N t (0) = 0)} d °- k 
k=l ^ ' 



oo do /j \ 

qa(2) = ^E^oEm 

ti =2 fc=l V 7 



E a ^ 

d fc =2 



[P t (N t (0) = 0)] 



do—k 



(2.29) 



(2.30) 



8 



In (gUID we already found that Sf k = -^e" 1 ^ - e~ (dfe ~ 1)4 ). So, we have 



do 



QA(2) = E^e- (rf0+1)t E(t 

d =2 k=l ^ 

oo do / , \ k 



X | JE. p -(d-l)t 



,d=2 



d=2 



d-l 



- E^E(t)E(J)(-D'(E^ 

d =2 fc=l v 7 i=0 v 7 \d=2 / 

/ oo ( d _ 1)t \ l 

= E °*E (t) E OV- 1 )' (Ed^i z < Vtf0+1) ' < 2 - 31 > 

d =2 k=l V 7 i=0 V 7 \d=2 / 

with Z t := o, c ie~^ d ~^ t /(d ~~ 1)- So, by integration we find 

QA(2) = f> E (t) E (") ("I)* (E " ' 1 Zle-^du (2.32) 



□ 



Example 

We would like to give an example where a closed-form solution is available which is free of 
integrals and gives us the time-dependent behavior of densities in the first and second line 
as sums whose main terms are exponentials in the time. Let us consider the special case 
where there are only two possible numbers of neighbors a and 6 on the random tree, i.e. 
we take G(s) = p a s a + p s b . We find 



Qpt(i) 



Pa 

a + 1 



1 - e 



6+1 



1 - e 



-(b+l)t 



(2.33) 



For the second layer we need to calculate the quantity Ct(n,x) := f Q Z™e ( x+1 ^ u du. We 
have 



Pa e -(q-l)t + Pb c -(b-l)t 

a — 1 6 — 1 

n\ ( Pa \ ( Pb \ [{a - 1){n ^ )+{h -l)j]t 



E 

3=0 



a - 1 



6-1 



(2.34) 



and so 



C t (n,x) 



_L_ (i _ e -(*+i)t) if n = 

<Ej=0 IjJ (a-l)(n-j)+(&-l)j+x-+l II > U 



(2.35) 



9 



So, for the second layers density we find the closed form 



Opt (2) = Pa 



1 



-(o+l)t fe-(a+l)t 



(a + 1) 2 (a + 1) 2 a + 1 

Pa 



fc=l v 7 i=0 



+ »eqe(: 

fc=l v 7 i=0 v 




6 + 1 



+ 1 



a- 1 6-1 



Pa Pb 



a-1 6-1 



k—i 



C t (i,b) 



(2.36) 



2.4 Procedure to Derive Higher Layer Densities 

It is natural to ask whether the procedure we just described to obtain densities on the first 
two layers can be generalized to obtain densities in a finite number of layers. To see the 
issue of higher layers more clearly let us specialize from the tree to the line. In this case, 
we claim that the time-dependent probabilities of the occurrence of any single-site pattern 
describing occupations up to a given finite height can in principle be calculated. However, 
in most cases the (probability of occurence of a) pattern can not be calculated directly but 
by a recursive algorithm which involves the computation of simpler patterns which we call 
the pre-image motives. 

The following procedure provides a method to find the time-dependent formula of the 
proportion of any pattern on a vertex. It consists of four steps: 1. find the pre-image mo- 
tives (the configurations from which the pattern under interest can increase or decrease), 
2. obtain the solutions of probabilities for occurence of the pre-image motives, 3. construct 
a differential equation of the target pattern based on the pre-image motives, and finally 4. 
solve the differential equation. 

As an example how the program works we will now calculate the probability of the occur- 
rence of Y t = (0, 1,0, l)' t , meaning the probability that the first and third layer are occupied 
and the second and fourth layer are empty. That the procedure stops after finitely many 
steps is not obvious from the beginning. Responsible for this fact is the screening. This 
will become clear in the example below. In a model without screening like [3] it is not true, 
and a corresponding recursion produces an infinite number of local motives. 



2.4.1 Step 1: Find the pre-image motives 

In this step we have to find the patterns whose occurrences contribute to an increase or 
decrease of our target pattern. In the case of Y t = (0, 1,0, 1)£ we find four pre-image 
motives, i.e. 







X 


A 






X 













X 











\,A 3 = 






1 








1 


h 





fx 


X 


A 





1 





1 










1 


<V 



fx x x\ 

1 

1 1 
\0 1 0/ 



(2.37) 
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This notation indicates on which position and layer a particle has been deposited (1) and 
where not (0), and where no particles have arrived so far (x). A particle is denoted by 
a 1 and empty positions that will remain empty due to blocking of neighbors or to the 
screening effect are indicated with a 0. Positions at and beyond which no particle has 
arrived so far are indicated with a x . Indeed, the proportion of the occurrence of Y t will 
increase with the proportion of both A\ and A 2 . In both patterns a particle is able to be 
deposited on the third layer and complete the pattern of Y t . The new particle can not be 
screened by particles in higher layers. On the other hand, the occurrence of A% or A^ may 
lead to a decrease of Y t , because they allow the arrival of a particle in the center location 
which results in (1,1,0,1)'. There are no other motives that can directly influence the 
proportion of Y t . 

2.4.2 Step 2: Obtain the solutions of the pre-image motives 

In this step we treat the pre-image motives one-by-one and find their solutions using the 
same four step procedure again. First we look at A\ and detect its pre-image motives. 

Finding A\ 

We apply the same procedure to find A\, With an abuse of notation we write A\(t) for 
the probabiity of its occurence. 

Step 1': The pre-image motives of A\ are: 

x x x x ' 

B x = | x x x x | ,B 2 = | x x x x | ,Ai = | 1 x| (2.38) 

10 10 

Step 2': Solutions of the pre-image motives of A\ 



B l {t)=te~ u (2.39) 








x 






i 









\o 


1 


0/ 



B 2 {t) = te- 6t P t (N t (l) = 0,m t (l) = l\N t (0) = 0) 



She~ 3t 



te~ 4t - te~ 5t (2.40) 



where we used our earlier result for Sf in (|2.17|) . 

Step 3': The differential equation for A\(t) takes the form 



A x {t) = B x (t) + B 2 (t) - 3Ax(t) 

= 2te- 4 *-te- 5 *-3A 1 (t) (2.41) 



Step 4': Solution of A^t) 
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Together with A^O) = we find 



Ai(t) = -e- 3t - 2e~ 4t - 2te~ u + -e" 5 * + -te~ 5t (2.42) 
4 4 2 

Finding A 2 

We apply the same steps to get A 2 . 

Step 1': The pre-image motives of A2 are 

(xxxx\ /xxxx\ /xxx\ 

x x 1 ,C 2 (t) = x x 1 ,A 2 (t) =1 1 (2.43) 
x 1 0/ \1 1 0/ \0 1 0/ 

Step 2': Now we solve Ci(t) and C 2 (i). 

(x x x x\ 
x x 1 = U t e- 2t . With f7 t := P t (N t (l) = 1, m t (l) = 1, N t {2) = 
x 1 0/ 
l,m t (2) = 2\N t (0) = 0). Now, we have to solve U t first. 

Ut = -2U t + Pt(N t (l) = l,N t (2) = 0,N t (3) = 0\N t (0) = 0) 

+P t (JV t (l) = l,iV t (2) = 0,iV t (3) = l,m t (3) = l\N t (0) = 0) 
= -2f7 t + te~ 3 * 

+ P t (iV t (l) = l,JV t (3) = l,m t (3) = l|JV t (0) = 0,JV t (2) = 0)e"* 
= -2U t + te~ 3t + P t (JV t (l) = l|JV t (0) = 0,N t (2) = 0) 

xP t (JV t (l) = l,m t (l) = l\N t (0) = 0)e"* 
= -2f7 t + te" 3 * + te _ *5 t 2 e _ * 
= -2^ + te- 3 * + te~ 2 *(e-*-e- 2< ) 

= -2?7 t + 2te" 3 ' - te" 4 ' (2.44) 

So, we have to solve U t +2U t = 2te~ 3t -te~ u . The homo geneous solution is Ut kom — 
and a particular solution is Ut, pa rt = —2te~ 3t — 2e~ 3t + \te~ 4i + \e" At . For t = we have 
Uq = 0. So, this gives the general solution U t = \e~ 2t - 2te~ 3 * - 2e~ 3t + \te~ u + \e~ u . 
Therefore, we have d{t) = U t e- 2t = if e~ 4 * - 2te~ 5t - 2e~ 5t + \te~ %t + \e~ %t . Now, we 
treat C 2 (t). We find 

(X x x x\ 
x x 1 = UtSfe- 1 (2.45) 
10 10/ 

So, 

C 2 (t) = (l^e" 2 ' - 2te- 3 * - 2e- 3t + l -te~ u + i e " 4 ^ (e"* - e^e"* (2.46) 
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Step 3': 

A 2 (t) = 2Ci(t) + 2C 2 (t)-3A 2 (t) 

= 7e~ u - H^e" 5 * - 8te- 5t + 5e" 6 * + 6te~ 6 * - ^e~ 7t - te~ 7t - 3A 2 (t) (2.47) 

Step 4': 

So A%{t) = Ce~ 3t and A p 2 {t) = -7e" 4t + f e" 5 *+4te- 5 *-| e - 6 '-2te" 6 * + ^e- 7 * + ite- 7 * 
and the general solution becomes 

A 2 (t) = ^e~ 3t - 7e~ 4t + ^e" 5 ' + 4te" 5t - -e" 6 * - 2te~ 6 ' + -e~ 7t + -te~ u (2.48) 
48 4 3 16 4 v ' 

Finding A3 and ^4 

The last two motives are much less complicated compared with the former two, so we can 
treat them together in one time. 

Step 1': The motives of ^3 are A\ and ^.3 itself, whereas the motives of A 4 are A 2 
and A 4 . 

Step 2': We already solved A\ and A 2 above. 

Step 3': The differential equations that we need to solve are 

A 3 (t) = A(t) - 3A 3 (t) and A 4 (t) = A 2 (t) - 3A 4 (t) (2.49) 

respectively. 

Step 4': The reader is invited to check that the solutions are 

A 3 (t) = -— e" 3 ' + -te~ 3t + 4e" 4t + 2te~ 4t - -e~ 5t - -te~ 5t and (2.50) 
4 4 4 4 

A 4 (t) = -^ e - 3t + ^te- 3t + 7e- 4t -fe- 5t -2te- 5t 

+ e- U + \te~ U ~ ^ ~ ±te~ 7t (2.51) 

2.4.3 Step 3: Construct a differential equation 

The differential equation for Y t is 

Y t = 2A 1 (t) + A 2 (t)-2A 3 (t)-A 4 (t) (2.52) 

Indeed, the appearance of Y t can increase by A\(t) and by its mirror motive. So, it counts 
two times. Also A 2 (t) increases the proportion of Yt but only one time, because its pattern 
is symmetric. Decrease of Y t occurs when a particle parks on top of A 3 (t) or A 4 (t) where 
the former counts two times because its mirror pattern has the same effect. 
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2.4.4 Step 4: Solve the differential equation 

After some calculations we find, using Y = 0: 



"3* + ^ te -3t + 7e -4t + 2te -4t _ 121-5* _ 3 St 

144 40 2 

784 6 112 16 [2 - b6) 

It should be clear that the probabilities of any other pattern can be computed in a 
similar way. 
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17 






+ 27 6 " 
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2.5 Comparison Results 

Let us now come back to the behavior on the first two layers and conclude the paper with 
a discussion of comparison statements for densities. 

Theorem 3 Consider a regular tree T^. The limiting density in the first layer is higher 
than in the second layer for all d>2. 



d 



Proof: Let us denote k (d) := {d/{d-l)) d {§ (d _$ +d+1 , so p£,(2) = ELo^ 1 )* 0k(d)~ 
jd+Yp ■ One verifies that 0k+i(d) / 0k(d) < 1 so that k i— > 0k(d) is decreasing. Therefore, 
making use of the alternating nature of the sum, we have the bound lim^oo p d (2) < 
0o(d) — 0i(d) + 02(d) — d/(d + l) 2 . The proof is then concluded by seeing that 0o(d) — 
(3\(d) + 02(d) < (2d + l)/(d + l) 2 for all d > 2. After some algebraic manipulations 
we find 0o(d) — 0i(d) + 02(d) = (d/(d — l)) d • we nave to check that 

(d/(d-l)) d ^^± < ^ or equivalents (l - ±) d > ^Ii){2d+1) • Developing the left term 

into a series and truncating it, we also find that (l — > 4^- — ^~g^~ 2 ^ = ~ Qdf^ ■ 
With equality only in the cases of d = 2 and d = 3. Furthermore, it is clear that 
2<yd 6tl^ +1 ^ ^ (3d-i)(2d+i) ^ or d > 2. So, finally, by checking the cases d = 2 and d = 3 
directly we conclude that the density of the second layer is strictly dominated by the first 
layer density for all d > 2. □ 

In [6] the issue of comparing the behavior of the process on a regular trees with that 
on a random tree having the same number of nearest neighbors on the average was raised, 
and a number of results were given. In our situation, we have the following. 

Theorem 4 Consider the random trees S and T with probability generating functions Gs 
and Gt respectively. If Gs(s) > Gt(s), for all < s < 1 then the first layer density of S 
exceeds the first layer density ofT for all t > 0. 

In particular, the first layer density of the regular tree dominates the first layer density 
of the regular tree for all t > if d' < d. 
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Proof: According to Theorem [2] the density of the first layer on a random tree S with prob- 
ability generating function G(s) = Y^=2 a n sn is given by Qpf(l) = Ylh=2 a k 1 ~ S k+i * ■ 
Define 7 (t) := Qpf (1) - QpJ(l) = ££L 2 (a fc ~ h) ( ) ■ We have 7 (0) = 0. In case 

of G s (s) > G T (s), the time derivative of 7 (i) becomes jfi(t) = ££L 2 ( a fc ~ ^)e" (fc+1)t = 
e- t (G s (e- t ) - Gr(e-*)) > for all t > 0. " " □ 

Theorem 5 Consider a regular tree with d neighbors for each vertex, with 

d € {2,3,4, . . . }. If random tree S has an average number of d vertex neighbors, then for 

any t > 0, the density of the first layer on S is higher than on T^. 

Proof: This follows from Jensen's inequality, if we can show that, for any fixed t, the 
function k i— > jdij(l — e~( fc+1 ^) = gt(k) has non-negative second derivative with respect to 

k. Indeed, a computation yields that = 2 ffff )f (e( fc+1 ) f -l-(fc+l)t- > 0. 

So, gt(k) is a convex function. □ 
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